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Six different types of cancer (i.e., breast, lung, colorectal, esophageal, gastric, and thyroid cancer) have high 
rates of incidence or mortality worldwide. It has been shown that activation of de novo lipogenesis is an early 
and common event in the cancer microenvironment. In this study, we performed lipid imaging and profiling 
for 134 tissue samples from six different types of cancer using matrix-assisted laser desorption/ 
ionization-Fourier transform ion cyclotron resonance mass spectrometry, with 2,5-dihydroxybenzoic acid 
and l,8-bis(dimethyl-amino)naphthalene as matrices in the positive and negative ion modes, respectively. 
Multivariate statistical analysis coupled with lipid distribution images revealed that significantly increased 
levels of monounsaturated fatty acids and monounsaturated phosphatidylcholines relative to 
polyunsaturated fatty acids and polyunsaturated phosphatidylcholines were observed in the cancer 
microenvironment compared with the adjacent normal tissue. The immunohistochemical assay indicated 
that fatty acid synthase, stearoyl-CoA desaturase-1, and choline kinase a were up-regulated in the cancer 
microenvironment compared with the adjacent normal tissue. Our findings suggest that de novo lipogenesis 
was activated in six types of cancer to promote a biosynthesis of lipids with monounsaturated acyl chains and 
to suppress a biosynthesis of polyunsaturated lipids in the cancer microenvironment. 

Cancer is a leading cause of mortality worldwide, accounting for around 13% of all deaths 1 . Therein, breast 
cancer in females and lung cancer in males are the frequently diagnosed cancers and the leading cause of 
cancer death; colorectal cancer is the third commonly diagnosed cancer in males and the second in 
females; gastric and esophageal cancer are followed by the high incidence and mortality rates 2 , and thyroid 
cancer incidence has also dramatically increased in the past 30 years 3 . All these cancers affect people's health 
seriously, and great efforts have been made to identify the mechanism of their development and progression, as 
well as diagnostic biomarkers. Lipids as a kind of important metabolites participate in the self assembly of 
phospholipids to form biological membranes and are involved in cell differentiation and signal transduction 
as second messengers 4 . Emerging evidence indicates that cancer and diverse human diseases are strongly corre- 
lated with changes in the levels of lipids 5,5 . 

Previous studies have noted that increased de novo lipogenesis is an early and common event in the develop- 
ment of cancer 7,8 , and cancer-related lipogenesis requires more lipid production to meet the requirements of rapid 
cell proliferation. During the above-mentioned process, a key lipogenic enzyme {i.e., fatty acid synthase, FASN) 
catalyzes the terminal steps in the de novo biosynthesis of saturated fatty acids (SFAs) from acetyl-CoA s . SFAs are 
partly converted into monounsaturated fatty acids (MUFAs) by stearoyl-CoA desaturase-1 (SCD1), and then 
generate monounsaturated phosphatidylcholines (MUPCs) mediated by a series of enzymes {e.g., choline kinase 
a (CKa) and lysophosphatidylcholine acyltransferase (LPCAT) 4 ). On the basis of previous data about breast, 
lung, and thyroid cancer, a significant up-regulation of MUPCs was observed in the cancer microenvironment. 
These findings indicate that cancer progression is associated with de novo lipogenesis 1012 . However, no invest- 
igation of in situ fatty acids reveals the role of de novo lipogenesis in cancer microenvironment. 

Matrix- assisted laser desorption/ionization-mass spectrometry imaging (MALDI-MSI) can simultaneously 
map multiple different lipids while preserving the morphological integrity of the analyzed tissue 13,14 . MALDI 
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coupled with ultrahigh resolution and ultrahigh mass accuracy 
Fourier transform ion cyclotron resonance mass spectrometry 
(FTICR MS) can greatly improve reliability and validity of MSI 15 . 
In addition, matrix is an important factor for MALDI-MSI, and the 
most commonly used matrix for detecting lipids in the positive ion 
mode is 2,5-dihydroxybenzoic acid (DHB). Previous studies have 
noted that l,8-bis(dimethyl-amino)naphthalene (DMAN) is an ideal 
matrix for studying low-molecular-weight metabolites in the nega- 
tive ion mode without any interference from the matrix, which could 
be used for the detection of fatty acids in tissue 1617 . Both matrices 
could be coated on tissue sections by a sublimation apparatus to 
generate high-purity and homogeneous matrix film, resulting in a 
significant improvements in lipid signals 18 . Furthermore, a strategy 
of merging both negative and positive ion MSI data from serial tissue 
sections was reported to provide an augmentation in the amount of 
useful information, and multivariate statistical analysis would screen 
biomarkers to distinguish cancer and normal regions of tissue 15 . All 
these MSI improvements will be of benefit for the study of changes in 
lipid distribution in tissue. 

In the present study, we investigated direct effects of de novo 
lipogenesis, especially the distribution of lipid species, on tissue sam- 
ples from 134 patients with six types of cancer (i.e., breast, lung, 
colorectal, esophageal, gastric, and thyroid cancer, Table 1). Both 
DHB and DMAN matrices were used to acquire mass spectra in 
the positive and negative ion modes, respectively, which contribute 
to an increase in the number of lipids detected. A combination of MSI 
and immunohistochemical assay was performed to reveal the link 
between changes in the levels of lipids and the expression of lipo- 
genic-related enzymes in cancer microenvironment. Our findings 
suggest that de novo lipogenesis has been activated in six different 
types of cancer and that increased activities of FASN, SCD1, or CKa 
may be associated with increased levels of lipids with monounsatu- 
rated acyl chains and decreased levels of polyunsaturated species. 

Results 

MSI of six different types of tissue samples. To investigate the effect 
of cancer-related lipogenesis in tissue, MSI was performed in this 
study for direct analysis of lipids in six different types of cancer (i.e., 
breast, colorectal, esophageal, lung, gastric, and thyroid cancer). 
DHB matrix was used to ionize phosphatidylcholines (PCs) and 
sphingomyelins (SMs) in the positive ion mode and DMAN was 
used to ionize fatty acids (FAs), phosphatidylethanolamines (PEs), 



and phosphatidylinositols (Pis) in the negative ion mode. In 
addition, a sublimation method for matrix deposition was adopted 
to improve repeatability and sensitivity of MSI, with the best signal 
intensity. For DMAN matrix, time-dependent changes in the 
intensities of metabolites in tissue under high vacuum condition 
were observed in Supplementary Fig. SI A and SIB. It was found 
that the absolute intensities of the metabolites almost remained 
unchanged in less than 25 min under high vacuum condition. The 
sublimation time and temperature for the different amount of the 
matrices were also systematically optimized to obtain a homoge- 
neous coating on the surface of ITO slides (see Supplementary Fig. 
SIC online). 

Representative mass spectra of two adjacent tissue sections from 
breast cancer tissue were acquired in the positive (left column, 
Fig. 1 A) and negative (right column, Fig. 1 A) ion modes, respectively. 
It should be noted that significant differences (e.g., PC(34:1) at m/z 
798.5426, PC(38:4) at m/z 832.5818, and C 18:1 at m/z 281.2484) were 
observed between the cancerous area and the adjacent normal area in 
the positive and negative ion modes, respectively. 

Representative ion images of MSI from six different types of can- 
cer in the positive and negative ion modes, with their individual 
corresponding H&E staining, are shown in Fig. IB. For the H&E 
staining, the cancerous areas and adjacent normal areas of six cancer 
tissues are circled by red dotted line and blue dotted line, respectively. 
The intensities of the cancer-related MUPCs (i.e., [PC (32:1) + K] + 
at m/z 770.5 103, [PC (34:1) + K] + at m/z 798.5426, and [PC (36:1) + 
K] + at m/z 826.5732) in the cancerous area were significantly 
increased compared with those in the adjacent normal area, whereas 
the intensities of polyunsaturated PC (PUPC, i.e., [PC(38:4) + K] + at 
m/z 832.5818) and phosphatide acid (PA, i.e., [PA(40:5) + K] + at 
m/z 789.4818) in the cancerous area were significantly decreased 
compared with those in the adjacent normal area. As shown in 
Fig. IB (right column), differentially expressed lipids obtained by 
MSI in the negative ion mode are FAs, PEs, and Pis. It should be 
noted that different kinds of FAs have shown different biological 
behaviors in the cancerous area. MUFAs (i.e., C 16:1 at m/z 
253.2174 and C 18:1 at m/z 281.2484) in the cancerous area presented 
high intensities compared with those in the adjacent normal area, 
whereas PUFAs, PE and PI with polyunsaturated acyl chains (i.e., 
C 2 2A at m/z 331.2645, [PE(38:4)-H] " at m/z 750.5441, and [PI(38:4)- 
H]~ at m/z 885.5499) in the cancerous area had a lower level than 
those in the adjacent normal area. However, some exceptions were 



Table 1 | Clinical information 


of tissue samples of 


six types of cancer 










Organ sites 


Breast 


Colon 


Esophagus 


Lung 


Stomach 


Thyroid 


Tissue No. 


50 


6 


15 


26 


19 


18 


Sex (F/M) 


50/0 


4/2 


3/12 


9/17 


14/5 


15/3 


Age (mean ± SD, range) 


44.7 ± 13.6 


63.7 ± 7.2 


63.0 ± 5.9 


58.0 ± 6.9 


58.5 ±11.6 


44.3 ± 10.5 


20-82 


55-75 


55-76 


40-68 


28-75 


31-56 


Cancer type 














Ductal carcinoma 


32 












Adenocarcinoma 


6 


4 




13 


16 




Squamous carcinoma 






14 


10 






Papillary carcinoma 












16 


Follicular carcinoma 












1 


NA 


12 


2 


1 


3 


3 


1 


Grading 














Gl 




1 


6 


6 


2 




G2 




1 


8 


12 


8 




G3 




2 




2 


6 




NA 


50 


2 


1 


6 


3 


18 


Tissue 














Cancerous 


33 


5 


12 


26 


18 


1 1 


Adjacent normal 


32 


5 


9 


4 


16 


12 


NA: no information on histological diagnosis or clinical stages. 
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Figure 1 | MSI of six different types of cancer. (A) Representative mass spectra from the cancerous area and the adjacent normal area of breast cancer 
tissue in the positive (left panel) and negative (right panel) ion modes. (B) Ion images of lipids from breast, colorectal, esophageal, lung, gastric, and 
thyroid cancer in the negative (right panel) and positive (middle panel) ion modes. Optical images of H&E stained sections adjacent to the section used for 
MSI are shown in left panel, and the cancerous area was circled with red dotted line and the adjacent normal area was circled with blue dotted line. 
MSI data were acquired with the spatial resolution of 200 |rm and 1000 um in the positive and negative ion modes, respectively. 



also observed for the alteration in the levels of these lipids, for 
example, the levels of C 2 oa anQ PI(38:4) were increased in thyroid 
cancer, whereas no difference in their levels was detected between the 
cancerous area and the adjacent normal area for esophageal and 
gastric cancer. Detailed information on the identification of these 
lipids of interest is shown in Supplementary Table SI. 



Lipid profiling of tissue samples. To further assess changes in the 
levels of the above-mentioned lipids, lipid profiling was performed 
on the basis of 134 tissue samples from six different types of cancer. 
Five mass spectra in the cancerous area and the adjacent normal area 
at five different positions, respectively, for each of tissue samples 
were acquired to screen cancer- related lipids and to assess tissue 
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heterogeneity. Finally, a total of 164 metabolites (or variables), 
corresponding to the m/z values after isotopic deconvolution and a 
combination of[M + H] + ,[M+Na] + , and [M + K] + as one m/z value, 
were enrolled in the final data set for multivariate statistical analysis. 

PLS-DA score plot revealed the obvious differences between the 
cancerous area and the adjacent normal area from the patients with 
breast, colorectal, esophageal, lung, gastric, and thyroid cancer, with 
the predicted residual sum of square (PRESS) of 0.4603, 0.5379, 
0.3757, 0.5069, 0.4125 and 0.4018, respectively (Fig. 2A). Y-scatter- 
plot of VIP values versus FC values (cancer vs normal) of the vari- 
ables was also used to select important variables. For breast cancer, as 
shown in Fig. 2B, the important variables with extremely large VIP 
and FC values can be easily found in the upper left or upper right 
hand quadrant. As exemplified by PC(32:1) and SM(22:0), both of 
FC values were 3.06 and —3.56, with VIP values of 1.55 and 1.75, 
respectively, indicating that both were significantly increased and 
decreased in the cancerous area compared with those in the adjacent 
normal area. All important variables with VIP >1 for six types of 
cancer were extracted on the basis of six Y-scatterplots (Fig. 2B). 
Finally, 18 common variables labeled in blue were extracted, which 
were strongly associated with six types of cancer, and they are C 16:1 , 
Ci8:i, C 20: i, PC(32:1), PC(34:1), PC(36:1), C 20:4 , C 22:4 , C 22:5 , 
PC(38:4), PC(38:6), PA(36:2), PA(38:3), PA(40:5), PE(38:4), 
PI(38:4), SM(22:0), and SM(24:1). The change trends of these 18 
lipids between the cancerous area (red) and the adjacent normal area 
(blue) are shown in Fig. 2C. It is worth noting that the levels of 
MUPCs (i.e., PC(32:1), PC(34:1), and PC(36:1)) and MUFAs (i.e., 
Ci6:i» C 18:1 , and C 20: i) in the cancerous area were significantly 
increased compared with those in the adjacent normal area 
(Wilcoxon-Mann- Whitney test, p <0.05), whereas polyunsaturated 
lipids (i.e., PC(38:4), PC(38:6), PA(38:3), PA(40:5), PE(38:4), and 
PI(38:4)) and PUFAs (i.e., C 20: 4, C 22;4 , and C 22:5 ) in the cancerous 
area were remarkably decreased compared with those in the adjacent 
normal area (Wilcoxon-Mann- Whitney test, p <0.05) except thyr- 
oid cancer. Detailed values of VIP, FC, and p are listed in 
Supplementary Table S2. It should be noted that elevated levels of 
some polyunsaturated lipids in thyroid cancerous area may be attrib- 
uted to different pathological phenotype compared with other 
cancer. Interestingly, remarkably decreased levels of some monoun- 
saturated and saturated SMs (i.e., SM(24:1) and SM(22:0)) were 
observed in the cancerous area compared with the normal area. 

Correlation of characteristic lipids in six different types of cancer. 

To further confirm the association of the levels of lipids with specific 
fatty acyl chain length and saturation with cancer, Spearman 
correlation analysis was performed to identify their relationships. 
As shown in Supplementary Fig. S2, for breast, colorectal, and 
esophageal cancer, MUFAs (i.e., C 16: i, Cig.i, and C 20: i) were 
strongly and positively correlated with MUPCs (i.e., PC(32:1), 
PC(34:1), and PC(36:1)), whereas both were negatively correlated 
with polyunsaturated lipids. PUFAs (i.e., C 20: 4, C 22;4 , and C 22;5 ) 
were strongly and positively correlated with polyunsaturated lipids 
(i.e., PC(38:4), PC(38:6), PA(38:3), PA(40:5), PE(38:4), and 
PI(38:4)); for lung and gastric cancer, the above-mentioned 
correlations became slightly lower; and for thyroid cancer, there 
were no obvious correlations among the above-mentioned lipids in 
the cancerous area. 

Association of the ratios of MUFAs/SFAs, MUPCs/SPCs, and the 
sum of FAs (or PCs) with cancer. It has been found that during de 
novo lipogenesis in vitro, SFAs are partly and endogenously 
converted into MUFAs by the rate-limiting enzyme of SCD1, and 
then the resulting products are converted by a series of enzymes to 
generate MUPCs. To confirm these phenomena in vivo, the ratios of 
MUFAs/SFAs (i.e., C 16:1 /C 16:0 and C 18:1 /C 18:0 ) and MUPCs/SPCs 
(i.e., PC(32:1)/PC(32:0), PC(34:1)/PC(34:0), and PC(36:1)/PC(36:0)), 
as well as the sum of FAs (i.e., Ci 6: i, Ci&o. Ci 8: i, and Ci 8: o) and PCs 



(i.e., PC(32:1), PC(32:0), PC(34:1), PC(34:0), PC(36:1), and PC(36:0)) 
in the cancerous area and the adjacent normal area for six types of 
cancer were investigated, respectively. As shown in Supplementary 
Fig. S3, most ratios of MUFAs/SFAs in the cancerous area were 
significantly higher than those in the adjacent normal area except 
the ratio of Ci 6: i/Ci 6: o in thyroid cancer. The ratio of PC(32:1)/ 
PC(32:0) in breast, colorectal, esophageal, lung, and gastric can- 
cerous area was significantly increased compared with that in their 
corresponding adjacent normal area, whereas this ratio in thyroid 
cancerous area was decreased. The ratio of PC(34:1)/PC(34:0) in 
lung or thyroid cancerous area was significantly increased or 
decreased compared with that in their individual corresponding adja- 
cent normal area, whereas no difference in this ratio was detected 
between the cancerous area and the adjacent normal area for breast, 
colorectal, esophageal, and gastric cancer. The ratio of PC(36:1)/ 
PC(36:0) in breast, esophageal, lung, gastric, and thyroid cancerous 
area was significantly increased, whereas no difference in this ratio 
was observed between colorectal cancerous area and adjacent normal 
area. In addition, it is worth noting that a common phenomenon for 
six types of cancer is that sums of FAs and PCs in the cancerous area 
were significantly increased compared with those in the adjacent 
normal area. These findings may imply that de novo fatty acid 
synthesis was increased in the cancerous area for six types of 
cancer, and cancer-related pathway converting MUFAs to MUPCs 
may was also activated during cancer development. 

Expression of FASN, SCD1, and CKa. We have performed immu- 
nohistochemical assays to examine the expression of FASN, SCD1, 
and CKa for a total of 134 tissue samples, including 105 cancerous 
regions and 62 normal regions. Immunohistochemical images from 
six representative tissue sections adjacent to the tissue sections 
analyzed by MSI are shown in Fig. 3A. As shown in Fig. 3B, it is 
found that three enzymes were significantly overexpressed in the 
cancerous areas except for CKa in esophageal and lung cancer and 
SCD1 in lung cancer. It is worth noting that immunohistochemical 
images of overexpression of FASN, SCD1, and CKa in the cancerous 
areas were strongly correlated with elevated levels of MUPCs and 
MUFAs in the cancerous areas detected by MS images (Fig. IB) or 
profiles (Fig. 2C). 

Distinguish of six different types of cancer. PLS-DA score plot 
further revealed the difference in the levels of these lipids in six 
types of cancer. As shown in Supplementary Fig. S4A, breast and 
thyroid cancer could be successfully differentiated from colorectal, 
esophageal, lung, and gastric cancer based on the first two principal 
components (tl and t2), suggesting that the expression of some lipids 
is significantly different between the former two cancers and the 
latter four cancers. To test the variable importance in discrimi- 
nating six types of cancer, the top 18 significantly changed lipids 
are shown in Supplementary Fig. S4B based on their VIP 
projection (VIP>1). As shown in the heatmap of VIP projection 
(right panel, see Supplementary Fig. S4B online), the lipids (i.e., 
PC(36:3), PC(32:1), and SM(34:1)) presented high level in 
colorectal, esophageal, and gastric cancer. For breast cancer, the 
levels of PC(38:3), PC(38:4), C 22:4 , C 20:2 , C 20:3 , PC(36:0), PC(38:5), 
and PC(36:1) were increased. The lipids of PC(38:2) and PC(34:1) in 
lung cancer and lipids of C 22:3 , PA(40:5), SM(36:2), PA(36:3), and 
PA(38:5) in thyroid cancer were also up-regulated. Taken together, 
our findings suggest that changes in the levels of the above- 
mentioned lipids may have cancer-specific metabolic features. 

Discussion 

Recently, emerging evidence indicates that the progression of cancer 
is usually associated with increased de novo lipogenesis 20 22 . In this 
study, we combined MSI and immunohistochemical assays to invest- 
igate the associations of cancer-related lipogenesis with the express- 
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Figure 2 | Metabolite profiling and feature selection of six types of cancer. (A) PLS-DA score plots of metabolites from the cancerous area (red dots) and 
the adjacent normal area (blue dots) of six types of cancer. (B) Y-scatterplots of VIP values versus FC values of metabolites from the cancerous area and the 
adjacent normal area of six types of cancer. 18 common variables labeled in blue have significant differences in their levels between the cancerous 
area and the adjacent normal area for six types of cancer. (C) Change trends of these 18 lipids in the cancerous area (red bar) and the adjacent normal area 
(blue bar) of six types of cancer (*, p <0.05; **, p <0.01). 
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Figure 3 | Immunohistochemical localization of the expression of FASN, SCD1, and CKa in six types of cancer tissues. (A) Representative 
irnmunohistochemical staining of FASN, SCD1, and CKa in six types of cancer tissues. Immunohistochemical assays of three enzyme were performed based on 
three serial tissue sections adjacent to their individual tissue section used for MSI and H&E staining in Figure 1, with two zoomed optical images (250X ) from the 
cancerous area and the adjacent normal area, respectively. (B) Scatterplots of the mean optical density of FASN, SCD1, and CKa expression. Quantitative analysis 
of the expression of three enzymes in the cancerous area and the adjacent normal area is shown as mean ± Standard deviation (*, p < 0.05; **, p < 0.01). 
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ion of lipids based on 134 clinical tissue samples from six different 
types of cancer. A MALDI-FTICR MS was employed to obtain the 
reliable distribution and accurate mass of lipids in clinical tissues. 
The use of both DHB and DMAN matrices provided an opportunity 
to obtain global analysis of lipid species on the basis of two adjacent 
tissue sections in both negative and positive ion modes, respectively. 
Based on the sublimation of DMAN, which had no matrix- related 
interference in negative ion mode, we have successfully performed 
fatty acid imaging or profiling for clinical tissues with well-preserved 
morphological features of tissues, suggesting that the different ion- 
ization modes can provide detailed information on metabolites in the 
cancer microenvironment. This is the first time that lipidomic 
changes (e.g., FAs, PCs, PAs, PEs, Pis, and SMs) have been system- 
atically characterised in conjunction with multiple lipogenic 
enzymes in six types of cancer tissue. It is worth noting that a con- 
sistent increase in the levels of the lipids with monounsaturated acyl 
chains where observed to be associated with the increased expression 
of FASN, SCD1, and CKa, and a relative decrease in the levels of 
polyunsaturated lipids in the cancerous area of clinical samples were 
observed. 

As shown in Fig. 1, MUPCs (i.e., PC(32:1), PC(34:1), and 
PC(36:1)) and MUFAs (i.e., C l6:l and Ci 8: i) were significantly up- 
regulated and species with polyunsaturated fatty acyl chains (i.e., 
PC(38:4), PA(40:5), PE(38:4), PI(38:4)) and PUFAs (i.e., C 20;4 and 
C 2 2-a) were down-regulated in the cancerous area of six different 
types of cancer compared with those in the adjacent normal area. 
These trends were further confirmed on the basis of multivariate 
statistical analysis (Fig. 2A) and Y-scatterplots (Fig. 2B) of lipid 
profile data from a total of 134 tissue samples. Finally, 18 common 
lipids (i.e., C 16:1 , C 18:1 , C 20: i, PC(32:1), PC(34:1), PC(36:1), C 20:4 , 
C 22:4 , C 22:5 , PC(38:4), PC(38:6), PA(36:2), PA(38:3), PA(40:5), 
PE(38:4), PI(38:4), SM(22:0), and SM(24:1)) with larger VIP and 
FC values were selected for six different types of cancer (Fig. 2C). 
These are consistent with previous studies about the overexpression 
of MUPCs in lung, breast, and colorectal cancer 10,12,23 . To further 
examine the relationship of these lipid abundances, Spearman cor- 
relation analysis of 22 lipids indicated a highly positive correlation 
between MUFAs and MUPCs and negative correlation between 
monosaturated and polyunsaturated lipids existed. These findings 
suggest that there may be similar lipogenic mechanisms to generate 
the above-mentioned lipids for six different types of cancer, but it 
should be noted that less correlation in thyroid cancer was observed, 
which may be associated with its lower aggressive behaviors com- 
pared with other cancers. 

Activation of de novo lipogenesis in cancer tissues that we focused 
on in this study was mediated by multiple lipogenic enzymes. FASN 
catalyzes the synthesis of palmitate (Ci 6: o) and stearic (Ci 8:0 ) acids 
from acetyl-CoA, with the overexpression in lung, breast, and pro- 
state cancer 8,24,25 , which is consist with our results. As shown in Fig. 3 
and Supplementary Fig. S2, increase in the sum of FAs (Ci 6: i, Ci 6: o. 
Ci8 : i, and C 18:0 ) was positively correlated with the up-regulation of 
FASN in six types of cancer, indicating that FAs synthesis may play 
an important role in cancer pathogenesis. In addition, increase in the 
ratios of MUFAs/SFAs (Ci 6: i/Ci 6:0 , C 18:1 /Ci 8:0 ) and MUPCs/SPCs 
(PC(32:1)/PC(32:0), PC(34:1)/PC(34:0), PC(36:1)/PC(36:0)) was 
also positively correlated with the overexpression of SCD1, which 
catalyzes the biosynthesis of MUFAs from SFAs, which are consist- 
ent with previous studies about lung, hypopharyngeal, prostate, 
ovarian, gastric, breast, and colorectal cancer 26 " 28 . 

Although the overexpression of CKa for the synthesis of PCs in the 
cancerous areas of breast, colorectal, gastric, and thyroid cancer was 
only observed, increase in the sum of PCs (PC(32:1), PC(32:0), 
PC(34:1), PC(34:0), PC(36:1), and PC(36:0)) was detected in all six 
types of cancer. These results may suggest that other enzymes maybe 
involved in the synthesis of PC, such as LPCAT4 23 , the expression of 
which was not determined due to sample limitations. In addition, 



PUFAs and polyunsaturated lipids were down-regulated in the can- 
cerous areas compared with the adjacent normal area. Taken 
together, the common lipogenesis pathway for six types of cancer 
was summarized in Fig. 4. Our results suggest that the use of in situ 
MSI detection coupled with immunohistochemical assay might pro- 
vide a new perspective to explore the common mechanisms of dif- 
ferent types of cancer. 

The development and progression of cancer can be viewed as a 
multistep process, with the complex interactions between deregu- 
lated lipogenesis and the cancer microenvironment 29 . It was found 
that increase in the degree of saturation could make cell membrane 
less fluid in most aggressive cancers 22,30 . Based on our observation 
that de novo lipogenesis would provide a consistent increase in the 
levels of lipids including monounsaturated acyl chains in the cancer- 
ous tissue, MUFAs and MUPCs replenish the membrane compo- 
nents for rapid cell proliferation and also pack the membrane more 
densely. This alteration of membrane fluidity may promote the car- 
cinogenesis in six types of cancer. On the other hand, polyunsat- 
urated acyl chains are more susceptible to peroxidation, and 
oxidized lipid species may be ultimately degraded into smaller react- 
ive products such as 4-hydroxyalkenal and malondialdehyde, which 
can cause cell damage 31,32 , which are consistent with our observation 
that modulation of de novo lipogenesis has increased the ratio of 
monounsaturated lipids to polyunsaturated lipids in the cancerous 
tissue compared with that in the adjacent normal tissue, cancer cells 
are protected from oxidative stress-induced cell death. More impor- 
tantly, different lipid metabolism may exist in thyroid cancer, result- 
ing in its less aggressive behavior. 

The development and progression for the six different types of 
cancer may involve different processes. As shown in Supplemen- 
tary Fig. S4, changes in the levels of 18 common lipids (i.e., 
PC(38:3), C 22:3 , PC(38:4), C 22:4 , C 20:3 , C 20:2 , PC(32:1), PA(40:5), 
SM(36:2), PC(36:3), SM(34:1), PC(38:2), PC(36:0), PC(38:5), 
PA(36:3), PA(38:5), PC(34:1), and PC(36:1)) for six different types 
of cancer were strongly correlated with the types of cancer. 
Polyunsaturated lipids in breast, colorectal, and thyroid cancer were 
significantly increased compared with those in esophageal, lung, and 
gastric cancer, while the levels of monosaturated lipids in esophageal 
and gastric cancers were increased compared with those in breast, 
colorectal, and thyroid cancer. For lung cancer, however, elevated 
levels of PC(34:1) and PC(38:2) were observed. These findings indi- 
cate that some different lipogenesis mechanisms among six different 
types of cancer may be involved. 

In summary, in the present study, based on the integrated strategy 
of MSI (or mass spectrometry profiling) in positive and negative ion 
modes using both DHB and DMAN matrices, respectively, we suc- 
cessfully detected several kinds of lipids (i.e., FAs, PCs, PAs, PEs, Pis, 
and SMs) and investigated direct effects of de novo lipogenesis in the 
six types of cancer. After merging ion images with immunohisto- 
chemical staining results, a significant change in the abundance of 
monosaturated lipids relative to polyunsaturated lipids and the over- 
expression of FASN, SCD1, and CKa were observed in the cancer 
microenvironment. The alteration in the degree of lipid unsaturation 
generated by de novo lipogenic enzymes in the cancer microenviron- 
ment may related to carcinogenesis. The improved MSI platform 
may be powerful tool to study cancer mechanisms and to screen 
novel diagnostic biomarkers. 

Methods 

Tissue samples. All specimens were remaining tissues after clinical pathological 
experiments, and informed consent was obtained from each patient. The approval 
was given by the Ethics Committee of Institute of Basic Medical Sciences, Chinese 
Academy of Medical Sciences, and all experiments were performed in accordance 
with relevant guidelines and regulations. In this study, a total of 134 primary tumor 
tissues were recruited, which were from 50 patients with breast cancer (invasive ductal 
and adenocarcinoma), 6 with colorectal cancer (adenocarcinoma), 15 with 
esophageal cancer (squamous carcinoma), 26 with lung cancer (squamous carcinoma 
and adenocarcinoma), 19 with gastric cancer (adenocarcinoma), and 18 with thyroid 
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Figure 4 | The lipid metabolic network of significantly changed lipids in six types of cancer tissues. The key enzyme of FASN catalyzes the synthesis 
of SFAs from acetyl-CoA. SFAs are endogenously converted into MUFAs through the part influence of SCD1, and then MUPCs are generated from 
MUFAs mediated by a series of enzymes (e.g., CKa). PUFAs could be produced by desaturation and chain elongation, followed by involving in the 
synthesis of polyunsaturated lipids (e.g., PUPCs). Red and blue arrows show up- and down-regulated expressions, respectively, in the cancerous areas of 
six types of cancer. Red and blue background indicates the up- and down- regulated levels of lipids, respectively. Ion images from breast cancer tissue are 
shown in the up panel. 



cancer {papillary and follicular). Cancerous area and adjacent normal area for each of 
tissue samples were carefully distinguished by the microscopic examination. The 
detailed information is listed in Table 1. 

Sample preparation and matrix deposition. All tissue samples were snap-frozen in 
liquid nitrogen and stored at — 80 "C until use. A 12 urn thin tissue section was cut and 
then thaw-mounted onto a glass slide coated by indium-tin-oxide (ITO) in a freezing 
microtome. After the resulting slides were rewarmed at room temperature, DHB or 
DMAN (Sigma-Aldrich) was coated using the sublimation apparatus {a homebuilt 
device). The sublimation was operated under the fixed vacuum (5 mTorr) at 120°C 
for 5 min. Finally, approximately 1 1 mg of DHB or 25 mg of DMAN was coated on 
each glass slide, and then recrystallization was performed to enhance signal intensity 
as described previously 33 . 



above-mentioned mass ranges. For the profiling data acquisition, the experimental 
parameters were the same as those of the MSI data acquisition. Five mass spectra at 
different positions of the cancerous area and the adjacent normal area were randomly 
acquired to confirm data reliability and to assess the heterogeneity of tissues. 

Structure identification of differentially expressed lipids of interest. Identification 
of lipid species of interest was first performed based on the Lipidmaps (http://www. 
lipidmaps.org/tools/index.html), with a mass tolerance of <2.5 ppm. Differentially 
expressed lipids were further identified on the basis of their measured accurate 
molecular masses, reliable isotope distributions, and tandem mass spectral data (see 
Supplementary Table Si online). The collision -induced dissociation experiments 
were selected to get tandem mass spectra of the lipids of interest. 



MSI or mass spectrometry profiling. MSI or mass spectrometry profiling of tissue 
sections was performed using an 9.4 T Apex-ultra rM hybrid Qh-FTICR MS (Bruker 
Daltonics, Billerica, MA, USA) equipped with a 355 nm Nd:YAG Smartbeam II laser 
(200 Hz, spot diameter 50 um). The experimental condition was optimized to obtain 
the best signal-to- noise ratio and good reproducibility in both positive and negative 
ion modes. For the MSI data acquisition, two adjacent tissue sections were analyzed at 
the spatial resolution of 200 um and 1000 urn in the positive and negative ion modes, 
respectively. A mass spectrum (or one pixel) was accumulated with three full scans 
with 50 and 70 laser shots each over a mass range of 400 — 1000 Da and 100 ~ 
1000 Da at the resolutions of 530,000 and 130,000 at m/z 400 in the positive and 
negative ion modes, respectively. As DMAN is easily sublimed under high vacuum 
conditions, low spatial resolution was used to ensure the reproducibility of MSI for an 
entire tissue section. Shown in Supplementary Fig. SI indicates that DMAN matrix 
can generate reproducible signals during the experimental time period of less than 
25 min. External calibration was performed using a lipid mixture (m/z 622.44423, 
m/z 734.56943, m/z 790.63203, and m/z 898.72593, which were from Avanti Polar 
Lipids, Inc., and identified by mass spectrometry) in positive ion mode and a tuning 
mixture (Part No. G2432A, Agilent Technology, Inc.) in negative ion mode over the 



Data handing and statistical analysis. Molecular images were visualized using the 
Flexlmaging software (version 2.1, Bruker Daltonics) with the absolute intensity. 
Peaks with signal-to-noise ratio of >3, relative intensity of >0.1%, and absolute 
intensity thresholds of 10,000 were selected as reliable variables, along with reliable 
isotope distributions using DataAnalysis 4.0 software (Bruker Daltonics). After 
isotopic deconvolution, the metabolites among different samples and different pixels 
were aligned within a narrow mass tolerance window (±0.001 Da) as one metabolite, 
and [M + H] + , [M+Na] + , and [M+K] + ions were also combined as one metabolite 
(or one variable). The intensities of all variables from one mass spectrum were 
normalized to a constant number of 1000, and the half of the baseline strength in each 
spectrum was adopted as their intensity of missing metabolites. Resulting datasets in 
the positive and negative ion modes from two adjacent tissue sections were combined 
into one dataset, and then transfer into Microsoft Excel for the following statistical 
analysis. SAS software (version 9.2, SAS Institute Inc., USA) was used for partial least 
squares discriminant analysis (PLS-DA), Wilcoxon- Mann- Whitney test, and 
Spearman correlation analysis. Metabolites with p value of less than 0.05 (Wilcoxon- 
Mann- Whitney test) were considered to be statistically significant and variable 
importance in the projection (VIP) of greater than 1.0 had differential reliability. 
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GraphPad software was used to generate box-chart. Important variables were further 
selected using a Y-scatterplot of VIP values versus fold change (FC) values. 

Immunohistochemical assay. Three lipogenic- related enzymes {i.e., FASN, SCD1, 
and CKa) were chosen for analysis by immunohistochemical staining on the basis of 
the fact that they were related to synthesis of SRAs, MUFAs, and MUPCs, respectively. 
Briefly, three serial frozen tissue sections, which were adjacent to the tissue sections 
analyzed by mass spectrometry, were warmed at room temperature for 30 min. After 
quenching of endogenous peroxidase activity with 3% H 2 0 2 for 10 min, theses 
sections were blocked with bovine serum albumin for 60 min at room temperature. 
These resulting sections were first incubated at room temperature with antibodies 
against FASN (1 : 250, Abeam, Cambridge, MA), SCD1 (1 : 200, Abeam), and CKa 
(1 : 200, Abeam), respectively, and then stored at 4°C overnight followed by 
rewarming at 37 C for 45 min. After washing with PBS three times, these sections 
were allowed to react with a secondary antibody against rabbit IgG (1 : 800), mouse 
IgG (1 : 200), mouse IgG (1 : 200), respectively. Detection was performed with the 
EnVision horseradish peroxidase system. To reduce nonspecific staining, phosphate 
buffer was used as a negative control. Adjacent normal area from tissue samples was 
used as an internal negative control and the cancerous area was used as an internal 
positive control for immunohistochemical staining. 

All immunohistochemical photographs were analyzed by using Image ProPlus 
(IPP, version 6.0, Media Cybernetics, Silver Spring, MD, USA) as described prev- 
iously 34 . In each area, five views were randomly selected under a microscope (250X). 
The mean optical density, as a quantitative measurement of stain intensity, was 
computed to determine average enzyme expression. 
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